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METHODS AND SYSTEMS FOR REDUCING DEPTH OF FIELD OF HYBRID 

IMAGING SYSTEMS 

STATEMENT REGARDING FEDERALLY SPONSORED RESEARCH 

[0001] The United States Government has certain rights in this invention pursuant to 
contract number DAAD 19-00-1-0514 awarded by the U.S. Army Research Laboratory in 
conjunction with the U.S. Army Research Office. 

RELATED APPLICATIONS 

[0002] This application claims priority to U.S. provisional application serial number 
60/417,332, filed 09 October 2003 and hereby incorporated by reference. 

FIELD OF THE INVENTION 

[0003] The invention relates generally to imaging systems and, more particularly, to 
hybrid imaging systems with reduced depth of field. 

BACKGROUND 

[0004] In prior art imaging systems, the image of an object gradually goes out of 
focus as the object moves from the plane of best focus. The image of a portion of the object 
that is in focus is also degraded by the image of adjacent portions of the object that are out of 
focus. This effect causes certain problems when detailed information is desired from the in- 
focus portion of the object and without influence from the surrounding portions of the object. 
It is particularly important in the field of microscopy to avoid degradation of the in- focus 
image due to adjacent out of focus images. 

[0005] Previous methods of obtaining a clear image of the desired portion or plane of 
an object include the use of pupil-plane filters. Pupil-plane filters utilize either amplitude 
(absorption) modulation or phase modulation of the light distribution in the pupil plane. 
Continuously varying amplitude pupil plates and annular binary pupil plates have been used 
to reduce the width of the central lobe of the axial intensity point spread function (PSF). 
These amplitude plate-based methods share two serious drawbacks: decreased optical power 
at the image plane and possible decrease in the lateral image resolution. A phase-only pupil 
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filter has also been used to reduce the axial spot size of a confocal scanning microscope. 
However, such a filter is not applicable to a hybrid imaging system because it employs a 
phase filter to reduce the width of the axial main lobe. However, due to the extremely high 
side-lobes in the PSF of such a phase filter, the useful optical power is reduced significantly. 

[0006] Structured illumination is another prior art approach to reducing the depth of 
field of an imaging system. For example, M. Neil et al., Method of obtaining optical 
sectioning by using structured light in a conventional microscope, Optics Letters, vol. 22, no. 
24, pp. 1905-1907 (1997), demonstrated that sinusoidal fringes of light would be formed by 
interference and projected onto the object. When an image is formed, the fringes go out of 
focus faster than a normal image. This effect leads to a slightly smaller depth of field; but the 
portions of the image that lie in the nulls of the sinusoidal fringes are lost. By the use of 
multiple exposures where the sinusoidal fringe is moved by a fraction of the fringe period for 
additional images, the complete image of the object can be retrieved when all of the images 
are superimposed on one another. One disadvantage of structured illumination is that precise 
alignment is needed. Another disadvantage is that multiple exposures must be made of the 
object to obtain a single image. This necessity for multiple exposures is problematic, 
especially when the object is moving, as in the case of live objects or moving parts along an 
assembly line. In the case of fluorescence microscopy, since the fluorophore is being 
bleached by the ultraviolet light used to excite the fluorophore, the later images are dimmer. 
If the object is a live cell, the ultraviolet light also damages the cell, making additional 
exposures particularly harmful. 

[0007] In confocal microscopy, optical "slices" are produced by focusing a point 
source onto the specimen and by imaging that point onto a point detector. Out-of-focus light 
is preferably removed to produce an in- focus image. Unfortunately, in order to obtain a 
complete image of the specimen, each plane of the specimen must be scanned point by point 
and the images of each plane then combined to achieve a three-dimensional result. 
Therefore, confocal microscopy is time consuming and is not suitable for imaging rapidly 
changing objects - such as living or moving samples. 

[0008] In contrast, the deconvolution microscope works with images of slices taken 
by a standard imaging system. Once images of the slices are taken and stored, along with the 
in-focus and out-of-focus PSFs of the imaging system, the deconvolution microscope 
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attempts to calculate the image in each plane. However, the digital post-processing used in 
calculating the image generates many undesirable artifacts. 

SUMMARY OF THE INVENTION 

[0009] The methods and system described herein support hybrid imaging with 
reduced depth of field, thereby reducing or eliminating the foregoing problems in the prior 
art. In one aspect, a three-dimensional hybrid imaging system is provided with increased 
optical-sectioning due to reduced depth of field and increased axial resolution, as compared 
to the prior art. Such an imaging system is, for example, useful in applications utilizing 
imaged planes of translucent three-dimensional objects. For example, the system has specific 
application with multi-level optical disk readers, optical medical (cellular) diagnostics, and 
microscopy. It is also useful for range determination. 

[0010] In one aspect, a hybrid imaging system is provided for imaging an object onto 
an image plane. A traditional imaging system has a defined value of axial resolution. The 
hybrid imaging system includes a modification of the traditional optical arrangement and is 
configured to form an intermediate image of the object at an intermediate image plane. A 
digital processing arrangement is configured to receive the intermediate image and to form a 
final image of the object. An optical element cooperates with the optical arrangement and the 
digital processing arrangement such that the axial resolution is increased over and above the 
defined value of the traditional imaging system. 

[0011] In another aspect, the optical element is configured to cooperate with the 
optical arrangement and the digital processing arrangement such that the depth of field is 
reduced below the defined value. 

[0012] In still another aspect, a method is provided for imaging an object onto an 
image plane, including the steps of: (1) forming an intermediate image of the object at an 
intermediate plane by one or more optical elements defining a first depth of field; (2) 
perturbing a phase function of electromagnetic energy forming the intermediate image by a 
specially designed optical element; (3) digitally processing the intermediate image to form a 
final image such that the one or more optical elements and the specially designed optical 
element define a second depth of field that is less than the first depth of field. 
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[0013] In one aspect, the electromagnetic energy comprises a spectrum of light, such 
as a visible light spectrum or an infrared light spectrum. The spectrum of light can have a 
range of incoherent wavelengths. 

[0014] In one aspect, the specially designed optical element comprises a rectangular 
phase plate with a two-dimensional plate function P(x,y). In one aspect, the phase function is 
randomized. In another aspect, the phase plate is a phase grating. In still another aspect, the 
specially designed optical element is formed by one or more optical elements and/or optical 
surfaces, in reflection, diffraction or refraction. 

[0015] In a further aspect, the method includes the step of configuring the specially 
designed optical element to cooperate with the step of digitally processing such that the final 
image formed from one portion of the object is in focus while the final image formed from 
adjacent portions of the object is out of focus. 

BRIEF DESCRIPTION OF THE DRAWINGS 

[0016] The present invention may be understood by reference to the following 
detailed description taken in conjunction with the drawings briefly described below. 

[0017] FIG. 1 is a schematic diagram of one hybrid imaging system. 

[0018] FIG. 2 shows a graph of an uncorrected zero-mean Gaussian random phase 
function with variance of p, to illustrate one phase function suitable for use within a phase 
plate in the hybrid imaging system of FIG. 1. 

[0019] FIG. 3 is a computer-generated graph of the magnitude of the ambiguity 
function of the random phase plate of FIG. 2. 

[0020] FIG. 4 is a computer-generated graph of the in-focus optical transfer function 
from the random phase plate of FIG. 2. 

[0021] FIG. 5 is an illustration comparing the best possible full-width half maximum 
of the hybrid imaging system of FIG. 1 to a standard optical imaging system in the prior art. 

[0022] FIG. 6 is a computer-generated profile of a rectangular reduced depth of field 
phase grating suitable for use in the hybrid imaging system of FIG. 1. 

[0023] FIGS. 7A-7D are computer generated graphs of the defocused, diffraction 
limited PSF of the hybrid imaging system of FIG. 1, with a rectangular reduced depth of field 
phase grating of FIG. 6. 
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[0024] FIGS. 8A-8D are computer-generated graphs of the defocused, diffraction- 
PSF of a standard imaging system in the prior art without the reduced depth of field phase 
grating. 

[0025] FIG. 9 is a computer-generated graph of the diffraction-limited Hilbert space 
angles for an optical imaging system with and without a reduced depth of field phase grating. 

[0026] FIGS. 10A-10D are computer generated graphs of the modulation transfer 
function (MTF) of a defocused, diffraction-limited imaging system that includes the reduced 
depth of field phase grating. 

[0027] FIGS. 1 1 A- 1 ID are computer generated graphs of the phase angles of the 
optical transfer function (OTF) of a defocused, diffraction-limited imaging system that 
includes the reduced depth of field phase grating. 

[0028] FIGS. 12A-12D are computer generated graphs of the MTF of a defocused, 
diffraction-limited imaging system without the reduced depth of field phase grating. 

[0029] FIGS. 13A-13H are computer generated images of a spoke target for a 
defocused, diffraction-limited hybrid imaging system with and without the reduced depth of 
field phase grating. 

[0030] FIG. 14 is a computer-generated graph comparing the diffraction-limited 
Hilbert space angles for a system with and without the reduced depth of field grating. 

r 

[0031] FIG. 15 is a computer- generated graph comparing the diffraction-limited and 
CCD-limited Hilbert space angles for an optical imaging system with the reduced depth of 
field phase grating. 

DETAILED DESCRIPTION OF THE PREFERRED EMBODIMENTS 

[0032] A hybrid imaging system combines a modified optical imaging system with 
digital processing; the modified optical imaging system is a standard optical imaging system 
modified with a specially designed optical element (hereinafter a "SDO element"). Unlike a 
standard imaging system that is cascaded with digital post-processing, in which the image 
acquired by the standard imaging system is manipulated in digital post-processing, the final 
image in the hybrid imaging system is obtained by digitally processing an intermediate 
optical image formed by the modified optical imaging system. Digital processing of the 
intermediate image makes the hybrid imaging system advantageous, among other reasons, 
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because additional degrees of freedom are provided to improve its imaging performance 
beyond the best feasible performance of a comparable standard imaging system. 

[0033] The SDO element is configured within the optical imaging system so as to 
produce an intermediate image which, with the digital processing, results in a final image that 
is sharp and in focus for a particular object range while other object ranges differ from the in- 
focus image more than in the standard optical imaging system. In one embodiment, the SDO 
element is a specially designed phase plate. In another embodiment, the SDO element is a 

combination of optical components designed to affect the phase of lightjransmitted 

■ 

therethrough in a specific way. The SDO element may be a refractive element, diffractive 
element, or a combination of refractive and/or diffractive elements. The SDO element may 
be a stand-alone physical element or one or more optical surfaces. The SDO element may be 
disposed in the optical imaging system at a location between the object being imaged and the 
final image plane; a detector typically resides at the final image plane to capture the image of 
the object. The SDO element may also be placed at or near to the aperture stop, or at or near 
to the image of the aperture stop, of the incoherent hybrid imaging system. 

[0034] FIG. 1 shows one hybrid imaging system 10. Hybrid imaging system 10 
includes a modified optical imaging system 12. Modified optical imaging system 12 has a 
lens 14 and a phase plate 16; phase plate 16 is positioned at the exit pupil of lens 14. 
Modified imaging system 12 accepts input light 20 from an object 1 1 and forms an output 
light 22 to form an intermediate image 24. Intermediate image 24 is then processed by a 
digital signal processor 30 to form a final image 34. 

[0035] Those skilled in the art appreciate that light 20 is electromagnetic radiation 
with a range of wavelengths, for example visible light (e.g., 400-750nm) or infrared light 
(e.g., 8-12 microns). Those skilled in the art also appreciate that phase plate 16 can be 
positioned at other locations, for example at the image of the exit pupil, at or near an aperture 
stop (not shown) of system 10, or at or near to the image (not shown) of the aperture stop. 
Further, those skilled in the art appreciate that lens 14 may be a single lens or represent a 
series of optical elements that cooperate to form image 24. 

[0036] Hybrid imaging system 10 does not substantially decrease optical power at 
the image plane by use of phase plate 16, as compared to use of an amplitude plate or a 
complex plate in the prior art. Phase plate 16 does attenuate the MTF of lens 14, thereby 
attenuating most frequency components in the image. Phase plate 16 may also distort the 
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phase of the OTF of lens 14. Accordingly, digital signal processor 30 in one embodiment 
incorporates a restoration digital filter that amplifies the attenuated frequency components of 
intermediate image 24 and, if necessary, corrects the phase of those frequency components. 
For example, the digital signal processor of one embodiment amplifies the attenuated 
frequency components of the intermediate image - and, if desired, corrects the phase of the 
components - with a linear frequency-domain inverse filter having a frequency response of: 



^Inverse {fx » fy ) ~ ' 



^Clear aperture {fx > f y) jj ( f /- \ , n 

— l f f \ > M Phase plate V x > J y ) * U 

" Phase plate \Jx>Jy) 



0 » ^Phase plate (/* > fy ) * 0 

0) 

where H C iear aperture (f X5 fy) is the in-focus OTF of lens 14 with a clear aperture (i.e., without 
phase plate 16 at its exit pupil) and Hp h asepiate (fx, f y ) is the in-focus OTF of modified optical 
imaging system 12 with phase plate 16 at the exit pupil. Since the inverse filter, H P h aS e plate (f x , 
fy), is a high-pass filter, it reduces the overall system signal-to-noise ratio. Those skilled in 
the art appreciate that nonlinear signal processing can be used to reduce noise. 

[0037] Although phase plate 16 is shown as a single element in FIG. 1, those skilled 
in the art appreciate that phase plate 16 may be a combination of optical elements designed to 
cooperate to yield the desired phase effect. Phase plate 16 may be refractive, diffractive or a 
combination of refractive and diffractive optical surfaces or physical elements. 

[0038] As described below, phase plate 16 can be designed in a variety of ways to 
reduce the depth of field of hybrid imaging system 10. In one embodiment, phase plate 16 is 
a rectangular phase plate with a two-dimensional plate function P(x,y). Since two- 
dimensions plate function P(x,y) is mathematically separable, between x and y, only one 
dimension is discussed below. 

RANDOM OPTICAL PLATE APPROACH 

[0039] The OTF of an incoherent optical system, as a function of defocus, can be 
shown as: > 

H{u,y,)= /^(p^ \u\ < 2 

where * denotes complex conjugate, and P(x) is the optical plate function: 
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P( X ) = 



e je(x \ for |*| <1 

0, otherwise 



(3) 

where J ~ and 

e{x) 

is some unspecified function. Since the magnitude of P(x) is unity, 
we are considering only phase plates that theoretically absorb no optical power. The defocus 
parameter \f/ is given by: 



nD 2 



^ = ^W 20 =kW 20 



1_ 1_ J. 

1/ d o d iJ 



(4) 

where D is the diameter of the lens and Xis the wavelength of the light. The distance do is 
measured between object 1 1 and the first principal plane of lens 14, while dj is the distance 
between the second principal plane of lens 14 and image plane 24. The quantity f is the focal 
length of lens 14. The wavenumber is given by k while the defocus aberration constant is 
given by W 2 o- 

[0040] The ambiguity function is given by: 

A(u, v) = $P(x + u 1 2)P* (x-ul 2)e J27rvx dx, \u\ < 2 

-(!-«/ 2) (5) 

The OTR and the ambiguity function are related by: 
H(u,y/)= A{u y uy/ 1 7t\ \u\ <r2 



(6) 

In other words, the projection of a radial slice of slope if//xof the ambiguity function onto the 
horizontal u-axis yields the OTF for defocus parameter \p. 

[0041] From Equations (5) and (3), we can consider the ambiguity function as a series 
of Fourier transforms of a function Qu(x) given by: 

Q u (x)=P(x + u/2)P*(x-u/2)=e J0{x \ |jc|<1, \u\ |*|<1,<2 

&(x) = e{x + u/2)-e{x-ul2\ \x\ < 1, w < 2 

It can be shown that, if the OTF is to change a maximum amount as a function of defocus 
parameter \f/, then the Fourier transform of Qu(x) may also change a maximum amount as a 
function of the Fourier variable. If Qu(x) is a stationary random variable, then the Fourier 
transform of Qu(x) may be uncorrected as a function of the Fourier variable. Equivalently, 
the transformed random variables (i.e., ambiguity functions) A(u,v) and A(u,v-h^) may be 



(7) 
(8) 
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uncorrelated for a given (u,v) and small y. Therefore, an uncorrelated ambiguity function 
leads to an incoherent optical system 10 with a reduced depth of field. 

[0042] If Qu(x) from Equation (7) is chosen as a continuous uncorrelated zero-mean 

Gaussian random variable with variance °" 2 , then the one-dimensional autocorrelation of the 
related ambiguity function is: 

(l-u/2) (l-u/2) 

EA(u,v)A*(u 9 v + y) = E J \Q U {x)Q* u (y)e i2 ™e- J2 " {v+ r )y dxdy, \u\ < 2 

-(l-«/2)-(l-u/2) ^ 

where E denotes expected value. Since Qu(x) is an uncorrelated random variable, then: 

a 2 for x = y 
0 otherwise 



EQu{x)Q:{y) = 



(10) 



(l-u/2) 

EA(u,v)A'(u,v + y) = a 2 ^e'^dx, \u\<2 

-(l-u/2) 



(11) 



• = a 2 (2 - u) sinc(2^(l - u 1 2)), \u\ < 2 

where sinc ( x ) = sin ( x )/ x . The one-dimensional autocorrelation of the random ambiguity 
function is then distributed with a sinc profile. 

[0043] The expected power of the random ambiguity function is given by: 

(l-u/2) (l-u/2) 

E\A(u,vf=E j jQ u (x)Q:(y)e j2mx e- j2 ^dxdy, \u\<2 

-(l-u/2)-(l-u/2) 



= <T 2 



(l-u/2) 

jdx, \u\ < 2 

-(l-u/2) 



= a 2 (2-u), \u\<2 (12) 

[0044] From the above arguments, it is shown that, if the function Qu(x) is a 
uncorrelated zero-mean Gaussian random variable, then the expected power of the 
corresponding random ambiguity function has a triangular profile in the u dimension and is 
constant in the v dimension. 

[0045] We can construct an approximately uncorrelated Gaussian random variable 
Qu(x) by selecting 0(x) from Equation (8) as a zero-mean Gaussian random variable with 
variance c? and autocorrelation r(r). From statistical theory, sums of Gaussian random 
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variables are Gaussian distributed. Therefore, from Equation (8) is also Gaussian. It 
can be shown that the autocorrelation of the transformed random variable eJ$ ^ is then: 
E e j9{x) e' j9{x ' r) = e~ 2 " 2(, - r(r)) (13) 

[0046] If the random variable ^Wis uncorrected, then , where is 

the Dirac delta function. The autocorrelation of e then becomes: 



E e JS{x) e' J9(x - r) =< 



1 for r = 0 

e~ 2a otherwise 



(14) 

Hence, when is uncorrelated and has a large variance; e jS ^ is also approximately 
uncorrected. From Equations (9) and (11), the resulting random ambiguity function will also 

be approximately uncorrelated, and the phase function will produce an incoherent 
optical system that maximally changes as a function of defocus parameter or object range. 

[0047] Turning now to FIG. 2, one random phase function suitable for use in 

phase plate 16 of hybrid imaging system 10 is shown. Phase function has been chosen 
to be an uncorrelated zero-mean Gaussian random phase function with variance of n. 

[0048] The ambiguity function corresponding to the phase function of FIG. 2 is 
shown in FIG. 3. As expected, the ambiguity function shown in FIG. 3 is uniformly spread in 
the v-axis direction while decreasing in the u-axis direction. This ambiguity function also 
appears very "noise-like." .Even to the eye, the correlation between adjacent points of this 
ambiguity function appears to be low. 

[0049] The in-focus OTF corresponding to the random phase plate of FIG. 2 is shown 
in FIG. 4. This OTF also appears to be noise-like, as desired. The level of this particular 

OTF is low, but it may be increased by forming the random phase function &[ x ) with a 

specific autocorrelation such that e jS ^ has a low-pass power spectrum. When is 

_/»9(jr) 

uncorrelated, e has an approximately white power spectrum. The digital processing used 
in processor 30 may then incorporate an inverse filter of this OTF, or another defocused OTF. 
In effect, system focus is altered not by movement of lens 14 but by digital filtering with the 
inverse of different defocused OTFs. 



10 



Matter No. 405834 



[0050] A succinct measure of the performance of a reduced depth of field hybrid 
imaging system 10 is the full width at half maximum (FWHM) amplitude of the resulting 
PSF. The FWHM measure can be applied to hybrid imaging systems when the PSF exhibits 
a main lobe. The FWHM measure of the reduced depth of field hybrid imaging system may 
sharply increase as a function of defocus parameter in comparison to a standard imaging 
system. An illustration of the best possible FWHM of the hybrid imaging system 10 of FIG. 
1 and that of a standard optical imaging system is shown in FIG.5. The width of the in-focus 
PSF of the standard optical imaging system has been normalized to unity. As expected from 
geometrical optics, the width of the PSF from the standard optical imaging system increases 
linearly beyond the value of its nominal depth of field. This value can theoretically be shown 

to approximately equal ^ = 1 ; this value is also very close to the simulated depth of field, 
using the FWHM criteria, of the standard optical imaging system. The width of the PSF from 
hybrid imaging system 1 0 sharply increases beyond a threshold defocus value. This 

threshold defocus value is a function of the variance and autocorrelation 

9{ X ) 

Hybrid 

imaging system 10 will therefore image at a plane 1 1 in object space with approximately 
diffraction limited resolution. Objects not at this plane will form maximally blurred images. 

PHASE GRATING APPROACH 

[0051] Referring briefly again to FIG. 1, hybrid imaging system 10, which includes 
phase plate 16 in modified optical imaging system 12, produces a blurred intermediate image 
24. A specific restoration digital filter may then be used to amplify the attenuated frequency 
components of intermediate image 24 - and, if desired, correct the phase of those frequency 
components - thereby obtaining an in-focus final image 34. The degree of digital restoration 
of an out-of-focus object depends on the similarity between the in-focus digital filter used and 
the out-of-focus digital filter required. 

[0052] The angle in Hilbert space between any two functions is a measure of the 
similarity between these functions. The smaller the angle between the two functions, the 
more similar the two functions. Therefore, rather than using the defocus parameter, this 
angle in Hilbert space between a defocused PSF and the in-focus PSF can be used as a metric 
to quantify the blurring of a defocused image. 



11 



Matter No. 405834 



[0053] Although Hilbert space is used in the following analysis, other spaces can also 

be used to attain similar results. This angle ® - ^ - n fe is defined for any defocus parameter 
value y as 



COS0 = 



where the inner product of the in- focus PSF 
defined as 



\h(«,0f 



and a defocused PSF 



(15) 



is 



(M".0)| 2 »W"»H 2 ) = JM">°)| 2 \h(u,y/fdu 

-00 

Also, the length in Hilbert space of the in- focus PSF is defined as 



(16) 



h(u,0f\\ = 



j\h(u,0f\h(u,0fdu 



(17) 



and the length in Hilbert space of the defocused PSF is defined as 



CO 

\\h{u,y/f\h(u,y/fdu 



(18) 



As a measure of the similarity between a defocused PSF and the in- focus PSF, the angle & is 
a more general metric to quantify the blurring of a defocused image than the defocus 
parameter V , which, for a given imaging system, is a measure of the defocus distance only. 

[0054] The PSF of a defocused imaging system with a rectangular aperture can be 
written as: 

2 



-1 -1 



ydxdy 



'-i J 



(19) 



where V * and ^ y are the defocus parameters in the directions x and y, respectively, *™ x and 

^max are tne half-widths of exit pupils in the directions x and y , respectively and k i s a 
constant. Since the defocused PSF of Eq. (19) is mathematically separable, we can restrict 
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our analysis to a one-dimensional defocused PSF. We introduce a phase plate /to (e.g., 
phase plate 16) at the exit pupil and we drop all multiplicative constants, yielding: 

\dx . 



jexpj^x 2 )-^ f{ x ) 



(20) 

[0055] For reduced depth of field, a phase plate function /to , which results in 
maximum image blurring at a slightly defocused plane (specified by relatively small defocus 

parameter value), is desired. To obtain the desired phase plate function /to for reducing 
the depth of field, Equations (15) and (20) can be combined and the optimization problem: 

(\h(u,0f ,\h(u, W ) x \ 2 ) 
min-jr n it 

' |A(«.or||A( U ,^r| (21) 

solved for a relatively small value of the defocus parameter V* . 

[0056] The optimization problem of Equation 21 can be solved for V* ~ * , for 

example, by assuming that the desired phase plate function /to is represented as a Fourier 

series. If the phase plate function /to is assumed to be represented as a Fourier series, then 
the PSF of the hybrid imaging system is consequently an array of narrowly-spaced spots. « 
Such an array of narrowly spaced spots would change its overall shape considerably as the, 
spots spread out and overlap due to propagation. If the location of the image plane is not to 
change, the phase plate does not have focusing power. For that case, we choose to represent 

the phase plate function /to by an odd Fourier series: 

sin(2^w vx) 

(22) 

Combining Equations (21) and (22), the optimum values of the fundamental spatial frequency 

v and the coefficients ^ n can be obtained numerically by using the method of steepest 
descent. 

[0057] The initial value of 0 in Equation (15) corresponding to a standard optical 

imaging system with a clear aperture and a defocus value of ^ x " 1 is 0.0601 radians. The 
value of ^ corresponding to a modified imaging system 12 with a sum of sinusoids according 
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to Equation (22) and with a number of coefficients N = 1 is 0.1284 radians. When the 
number of sinusoids is increased to N = 5 , the value of # increases to 0.1355 radians, which 
is a change of only 5.52% from the value of 0 corresponding to N = 1 . Therefore, for 
example, we restrict the number of phase grating coefficients to N = 5 . 

[0058] The values of the fundamental frequency v and the coefficients ^ n 
corresponding to an F/4 hybrid imaging system are shown in Table 1. Due to the previously 
discussed mathematical separability of the defocused PSF, the desired two dimensional, 
rectangular phase plate function for reducing the depth of field can be expressed as 

5 5 

/(*» y) = £ b n sin(2^i vx) + ]T b n sin (ijm yy) 

1 1 (23) 

We refer to the phase plate corresponding to f( x >y) 9 whose coefficients are shown in Table 
1, as the rectangular reduced depth of field (RDF) phase grating. The profile of the 
rectangular RDF phase grating is shown in FIG. 6. 



Table 1 . Rectangular RDF phase grating parameters 



v (cycles/mm) 


b x (//m) 


b 2 (^m) 


b 3 (//m) 


b 4 (//m) 


b 5 (^ra) 


1.0400 


1.1705 


-0.0437 


0.0271 


-0.0325 


-0.0700 



> [0059] FIGS. 7A-7D show the defocused PSF of an F/4 diffraction-limited modified 

imaging system 12 with a rectangular RDF phase grating 16, with parameters as shown in 
Table 1, at the exit pupil of system 10. 



Table 1 . Rectangular RDF phase grating optimum parameters 



v (cycles/min) 


b, 


b 2 


b 3 


b 4 


b 5 


1.0400 


1.1705 


-0.0437 


0.0271 


-0.0325 


-0.0700 



It may be noted that the variation with defocus in the shape of the PSF shown in FIGS. 7A- 
7D is greater than the variation with defocus in the shape of the PSF of a similar standard 
diffraction-limited system (without the rectangular RDF phase grating - i.e., "clear 
aperture"), as shown in FIGS. 8A-8D. 
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[0060] The angles in Hilbert space between the in- focus PSF and defocused PSFs of 
an F/4 diffraction-limited standard imaging system and a modified imaging system including 
the rectangular RDF phase grating are shown in FIG. 9. The parameters of the modified 
imaging system are shown in Table 1. It can be seen in FIG. 9 that, for all of the defocus 
parameter values in the range shown, the angle in Hilbert space between the in-focus PSF and 
defocus PSF of the modified system 12 with the rectangular RDF phase grating at its exit 
pupil has a greater value than the corresponding angle in Hilbert space for the standard 
imaging system. As a result, the shape of the PSF of a diffraction-limited modified imaging 
system including the rectangular RDF phase grating at its exit pupil shows greater variation 
with defocus than the shape of the PSF of a similar standard diffraction-limited imaging 
system. Furthermore, for lower defocus parameter values, the Hilbert space angle shown in 
FIG. 9 of the modified system is larger than the corresponding Hilbert space angle of the 
standard system by more than a factor of two. Therefore, a diffraction-limited modified 
system with the rectangular RDF phase grating at its. exit pupil has less than half the depth of 
field of a similar diffraction-limited standard system. 

[0061] Turning now to FIGS. 10A-12D, the MTF and OTF of the defocused F/4 
diffraction-limited modified imaging system with the rectangular RDF are shown in FIGS. 
10A-10D and in FIGS. 1 1A-1 ID, respectively, and the defocused OTF graphs of a similar 
diffraction-limited standard system are shown in FIGS. 12A-12D. A comparison of FIGS. 
1 1 A-l ID and FIGS. 12A-12D reveals that there is a greater variation with defocus in the 
phase of the OTF of the modified system than in that of the standard system. 

[0062] FIGS. 13A-13H show computer-simulated images of a spoke target for 

# 

different defocus values. The images of FIGS. 13 A, 13C, 13E and 13G show the spoke target 
images obtained using an F/4 incoherent diffraction-limited standard imaging system with 
different defocus values. The FIGS. 13B, 13D, 13F and 13H illustrate the computer- 
simulated images of the spoke target, for corresponding defocus values, obtained using an F/4 
incoherent diffraction-limited imaging system that was modified by including the rectangular 
RDF phase grating. In each column, the value of the defocus parameter varies from top to 
bottom from 0.0 (i.e., in focus) to 3.0. It can be seen in FIGS. 13A-H that the depth of field 
of the imaging system that was modified with the rectangular RDF phase grating is reduced 
in comparison to the depth of field of the standard system. That is, the image formed with the 
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modified imaging system changes more quickly with increased defocus parameter compared 
with the blurring of the image formed with the standard imaging system. 

[0063] Similar studies can be performed to examine the effect on the reduction of 
depth of field of using a CCD optical detector rather than an ideal optical detector. FIG. 14 
shows the angles in Hilbert space between the in-focus PSF and defocused PSFs of an F/4 
CCD-limited standard imaging system and a modified imaging system including the 
rectangular RDF phase grating. The parameters of the rectangular RDF phase grating used in 
this model have been optimized for use with the CCD-limited system. It can be seen in FIG. 
14 that, for lower defocus parameter values, the angle in Hilbert space between the in-focus 
PSF and defocused PSF of a modified system with a rectangular RDF phase grating at is 
pupil is larger than the corresponding angle of a standard system by more than a factor of 
two. Therefore, a CCD-limited modulated system including a rectangular RDF phase grating 
at its exit pupil has less than half the depth of field of a similar CCD-limited standard 
imaging system. 

[0064] FIG. 15 illustrates the angles in Hilbert space between the in-focus PSF and 
defocused PSF of an F/4 diffraction-limited imaging system modified with a rectangular RDF 
phase grating having the parameters shown in Table 1 . Also shown in FIG. 1 5 are the angles 
in Hilbert space for a CCD-limited modified system, including a rectangular RDF phase 
grating with its parameters optimized for the CCD-limited system. For the defocus parameter 
values in the range shown in FIG. 15, the Hilbert space angles obtained from the CCD- 
limited system are smaller than the corresponding angles obtained from a diffraction-limited 
system. Therefore, as expected, the variation with defocus in the shape of the PSF of the 
CCD-limited system is greater than the variation with defocus of a similar system without the 
phase plate. 

[0065] Hybrid imaging system 10 has advantages because, among other reasons, the 
depth of field is decreased. Therefore, images of slices of three-dimensional objects can be 
improved by using hybrid imaging system 10. Phase filter 16 can utilize for example the 
afore-described random optical mask or the rectangular RDF phase grating. Hybrid imaging 
system 10 is more efficient than, for example, a confocal microscope because, unlike the 
confocal microscope which must scan the specimen using a point source and then image each 
point onto a point detector, system 10 is capable of rapidly capturing each "slice" so as to 
quickly accumulate a three-dimensional final image with high axial resolution. Also, system 
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10 is faster than an imaging system with structured illumination, which requires more than 
one exposure for every image. Furthermore, since system 10 includes digital processing in 
the image acquisition procedure, rather than for manipulation of images acquired by an 
optical-only imaging system, system 10 provides superior final images with fewer artifacts 
compared to, for example, a deconvolution microscope. Hybrid imaging system 10 further 
facilitates imaging of thinner optical "slices" of the specimen with high lateral resolution as 
well as high axial resolution, resulting in a final three-dimensional image 34 with higher 
resolution. 

[0066] Since certain changes may be made in the above methods and systems without 
departing from the scope hereof, it is intended that all matter contained in the above 
description or shown in the accompanying drawing be interpreted as illustrative and not in a 
limiting sense. For example, another type of high pass-filter, such as a Weiner filter, can be 
added to system 10 to restore the intermediate image. Nonlinear filtering can be used. It is 
also to be understood that the following claims are to cover all generic and specific features 
described herein, and all statements of the scope which, as a matter of language, might be 
said to fall there between. 
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